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ABSTRACT 

Microscopic diffusion processes (such as radiative levitation and gravitational 
settling/thermal diffusion) in the outer layers of stars are important because they 
may give rise to surface abundance anomalies. 

Here, we compare radiative accelerations (grad) derived from the new Opacity 
Project (OP) data ^ with those computed from OPAL and some previous data 
from OP. 

For the case where we have full data from OPAL (carbon, 5 points in the p—T 
plane), the differences in the Rosseland mean opacities between OPAL and the 
new OP data are within 12% and less than 30% between new OP and previous 
OP data (OPl) ^. The radiative accelerations grad differ at up to the 17% level 
when compared to OPAL and up to the 38% level when compared to OPl. 

The comparison with OPl on a larger (p — T) space gives a difference of up to 
40% for grad{C). And it increases for heavier elements. The differences increase 
for heavier elements reaching 60% for Si and 65% for S and Fe. 

We also constructed four representative stellar models in order to compare 
the new OP accelerations with prior published results that used OPAL data. 
The Rosseland means overall agree better than 10% for all of our cases. For the 
accelerations, the comparisons with published values yield larger differences in 
general. The published OPAL accelerations for carbon are even larger relative to 
OP compare to what our direct comparisons would indicate. Potential reasons 
for this puzzling behavior are discussed. 



^These data will be made generally available via a server or on CDROM. All request should be sent to 
Claude Zeippen (Claude.Zeippen@obspm.fr), or to Anil Pradhan (pradhan@astronomy.oliio-state.edu). 

^The OPl data are the data present in the database at Centre de Donnees de Strasbourg (CDS, 
http : / /cdsweb .u-strasbg.fr/). 
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In light of the significant differences in the inferred acceleration rates, theo- 
retical errors should be taken into account when comparing models with observa- 
tions. The implications for stellar evolution are briefly discussed. The sensitivity 
of grad to the atomic physics may provide a useful test of different opacity sources. 

Subject headings: Atomic data - opacity - radiative accelerations - diffusions - 
stars: interior - stars: evolution 



1. Introduction 

Element segregation processes arc clearly seen in stars and solar models and occur from 
well understood physical processes. Gravitational settling and thermal diffusion will tend 
to make heavier species sink relative to the light ones. Radiative pressure will tend to 
cause some species to rise in the stellar interior (we refer to this as radiative acceleration). 
Radiative accelerations have been calculated using the Opacity Project (OP, Seaton et al. 
1994, The Opacity Project Team 1995, 1997) data (Alecian & Artru 1990; Alecian, Michaud 
& Tully 1993; Alecian 1994; Gonzalez et al. 1995; Leblanc & Michaud 1995; Hui Bon-Hoa et 
al. 1996; Seaton 1997,1999, Alecian & Leblanc 2000; Leblanc & Alecian 2004). Evolutionary 
calculations, however, have almost all been based on the OPAL theoretical opacities (Iglesias 
& Rogers 1996) except for Seaton (1999). 

In this paper we compare results from the up-dated data from opacity project (OP) 
with those obtained from OPAL and previous data from OP. We begin by discussing the 
astrophysical impact of element separation processes, and then move to our motivation for 
comparing with other datasets. 

It is now generally accepted, as proposed by Michaud (1970) for Ap stars and Watson 
(1970) for AmFm stars, that radiative levitation plays an important role in hot and slow- 
rotating stars. The morphology of the Horizontal Branch of different globular clusters (GCs) 
presents some features that arc not predicted by standard stellar models: Gaps in the blue 
tail (Ferraro et al. 1998), jumps in the Stromgren color-magnitude diagram (Grundahl et 
al. 1999), surface gravity anomalies (Moehler et al. 1995) and abundance anomalies (Behr 
et al. 1999; Behr et al. 2000a,b). This results in a bi-modal distribution in the HB stars 
of the studied GCs. Qualitatively, Hui-Bon-Hoa and co-authors (2000) showed that these 
observations could be the signatures of radiative acceleration (grad)- 

Microscopic diffusion can also affect the internal structure of stars. The impact of 
gravitation setthng on solar models has been extensively explored (Bahcall & Pinsonneault 
1992,1995; Vauclair 1998, Turck-Chieze et al. 1998, Sackmann & Boothroyrd 2003). Diffu- 
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sion deepens the solar surface convection zone, improving the agreement with hehoseismic 
data on its depth, and yields a surface helium abundance in good agreement with the value 
deduced from hehoseismic studies. Element separation processes will affect the thermal struc- 
ture, the convection depth, and the inferred initial abundances of solar models. (Charpinet 
et al. 1997; Turcottc ct al. 1998, hereafter T98a; Turcotte, Richer & Michaud 1998, hereafter 
T98b; Richer et al. 1998, hereafter R98). 

In addition, globular cluster age estimates are affected by the inclusion of settling and 
diffusion (Chaboyer et al. 1992; Chaboyer, Sarajedini & Dcmarquc 1992; Vandcnbcrg et 
al. 2002) leading to an age reduction of order 10% relative to models that neglect settling. 
Michaud et al. (2004) showed the importance of micro-diffusion in the age determination of 
open clusters and its effect on isochrone morphology. It has become clear that radiative levi- 
tation and diffusion processes must be included in stellar evolution codes in a self-consistent 
manner. 

Computing these effects has been challenging; detailed comparisons between theory and 
observations have had mixed success. T98b found that their predicted overabundances were 
larger than the observed ones. They explored other physical processes, such as turbulent 
diffusion and mass loss (T98b; Richer et al. 2000), as potential solutions. However, it is also 
possible that the uncertainties in the diffusion velocities themselves could be a significant 
error source. Because the two effects (gravitational settling and Icvitation) are in opposite 
senses, a small difference in the radiative accelerations could change the magnitude (or even 
the sign) of the predicted abundance anomalies by a much larger than linear factor when 
the two effects are of similar magnitude. Before any attempt to generate a model with 
micro-diffusion including radiative levitation, a careful study of the atomic data available is 
necessary. We therefore evaluate the uncertainties in the radiative acceleration coefficients 
themselves. 

In section 2 we describe the method used to obtain the radiative accelerations and the 
different stellar model structures. In the first part of section 3 we compare and analyze the 
results from the new data and from the previous OP data. In a second part of this section 
we compare the results from OPAL and NEW for carbon. This is the only element for which 
we have some monochromatic opacities from OPAL. For the other elements we compare 
our results with published works in the third part of this section. The discussion and the 
conclusion constitute sections 4 & 5 respectively. 
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2. Method. 

In a star, the net radiative flux is going outward. As the photons move toward the 
surface, they interact with the ions present in the star. Each species, depending on its state 
of ionization, the excitation level and its corresponding absorption cross section, experiences 
a net force that tend to make ions rise when they absorb momentum from photons. The 
resulting acceleration is called the radiative acceleration and is defined as follows (Seaton 
1997): 

9rad{k) = -J^'^R^ik) (1) 

where M is the mean mass per atom, M(k) is the atomic mass of atom k, kr is the 
Rosseland mean, F is the total flux of the radiative source (a blackbody at T = T'e//) and 
^ corresponds to the total momentum radiative flux associated to it. Finally, 7(A;) is a 
dimensionless quantity characterizing the individual contribution of element k to the total 
opacity defined as follow (Seaton 1997): 

Ki,{t0tal) = ^ K^{k) + Kscat (3) 
k 

where K,^{k) is the monochromatic opacity of element k, K^{total) is the total opacity, 
both in cm^/g of mixture, and f,^ is a weighting function defined below. As a reminder, with 
this notation, kr is defined as: 



1 



Kr J Ky{total) 



Udv (4) 



■'^ ~ An^klT^ (eWfcB _ 1)2 

In the present work, we compute g^ad with and without including the effect of momentum 
tranfer to electrons during photoionization processes. In other words, we have use Ki,{k) or 
K^^°'{k) (as define in Seaton 1997, equation (31)). This allow us to isolate each difference in 
the atomic data. For the OPAL data as for the OP data we subtracted the electron scattering 
opacities and the opacities corresponding to the momentum transfered to the electrons during 
the photoionization from the monochromatic opacities of the studied elements. However This 
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is not subtracted from the total opacity. We used exclusively the OP data to remove these 
contribution. We did not remove the OPAL electron scattering opacity from the OPAL data 
because it is already subtracted via the OP data. 

In order to estimate the systematic errors in the radiative accelerations due to atomic 
data, we compare the values of grad obtained with four sets of atomic data, two from OPAL 
and two from OP (NEW and OPl). Wc could compare with OPAL data for carbon (Iglesias 
2004, private communication) at the five points in the (p-T) plane listed in Table 1, and 
discuss indirect comparisons with the literature in §3. Ideally, one would directly compare 
the monochromatic opacities (total and specific for each clement) from OP and OPAL. 
However, due to technical problems the OPAL group could not provide a full set of data. 

We first calculated the Rosseland mean opacities and 7(C) with the 3 sets of data. We 
then calculated the accelerations for different sets of Tgff and radius (these determine the 
flux in eq.l). 

For the other elements present in the mixture, we do not have the relevant OPAL atomic 
data. In section 3.2 we thus compare OP data for other elements with previously published 
OPAL results (T98a; T98b and R98). We calculated the Roseland mean opacities, 'j{k) 
and the radiative accelerations using OP data for different types of stars. This indirect 
comparison allows us to span a large portion of the (p-T) plane. As we will show, however, 
it is difficult to directly interpret differences obtained in this manner. 

We calculate the structure of different stars with the YREC code (see Bahcall, Pin- 
sonneault & Basu 2001 for a description of the model ingredients) in which we modifled 
the micro-diffusion subroutines. We extended the composition vector to include all species 
present in the mixture in order to track individually their evolution within the model as a 
function of time. The diffusion equation for element separation processes has been modi- 
fied in order to treat separately each species of the initial mixture. Instead of treating the 
gravitational accelerations of all element as if they were settling like fully ionized Fe, the 
gravitational settling coefficient are calculated individually for each element. We also in- 
cluded the individual radiative accelerations calculated using the OP data and the method 
described by Seaton (1997, 1999). All the technical details of the calculation will be presented 
in a subsequent paper. 

We first calibrated our model with a mass of 1.0 Mq to reproduce the observed solar 
luminosity (Lq), radius (Rq) and surface Z/X ratio at the solar age (4.57 Gyr). The initial 
composition used here includes 17 elements (H, He, C, N, O, Ne, Na, Mg, Al, Si, S, Ar, Ca, 
Cr, Mn, Fe, Ni) with relative abundances from the Grevesse & Noel (1993) mixture. The 
cahbrated model yields the initial helium mass fraction (Y© = 0.271) and the mixing length 
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(cKo = 2.09). We used this solar calibration for our other models. As previously noted (T98a), 
it might not be justified to use these parameters from a calibrated Sun for stars with a range 
of masses and evolutionary states. But by doing so we transform the absolute errors in the 
input physics into relative ones, which is a significant advantage for a solar calibration. For 
all models as well as for the 5 special points for which we have the monochromatic opacities 
from OPAL, we used the composition given in Table 2. 

In the present work we focus primarily on quantifying the errors in the coefficients 
themselves, so we have not yet included the feedback on the opacity induced by changes in the 
heavy element mixture. Changes in the relative heavy element abundances can certainly be 
important in some contexts as suggested by Alecian, Michaud & TuUy (1993). For example. 
Richer, Michaud & Turcotte (2000) found an iron convection zone in AmFm models including 
micro-diffusion that was caused by feedback effects. We intend to include such effects in 
future work focused on the applications of diffusion to stellar evolution problems. 



3. Results 

Our aim in this section is twofold: we present the effect of the up-dated OP opacities on 
the radiative acceleration and then compare the new results to OPAL. We present two specific 
cases, C and Fe, to show how the monochromatic opacities modify the radiative acceleration. 
The improvement in the OP atomic data due to the inclusion of inner-shell transitions 
systematically enhance the radiative accelerations relative to the older OP data. The impact 
of the new physics increases for heavier species and higher temperature. When directly 
compared to the OPAL data, the accelerations for carbon are in reasonable agreement. 
However this agreement deteriorate in the comparison with published work. The agreement 
between OP and OPAL is expected to be, and is, less favorable for heavier elements. 



3.1. Comparison between New OP and OPl. 

In order to understand the contribution of different ingredients to Qrad, we present the 
monochromatic opacities for C and Fe at a specific {p—T) point. The main difference between 
the two sets of data (OPl and New OP) is the inclusion of the inner-shell transitions in the 
later. 

For our purpose, it is worth recalling that the acceleration Qj-ad oc x 7 and 7 oc 
J i^eianent _ ^jj-^g lutcrplay between the various terms in equation (1) is illustrated in Figure 
1. The data used for this figure are from the new OP data and OP data without inner- 
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shell transitions (Seaton, private comunication) . The top panels show the monochromatic 
opacities for C (left) and Fe (right). The middle compare the total monochromatic opacities 
(same on each side) and the bottom panels illustrate the ratio of the 2 which corresponds to 
part of the integrand in the definition of 7 (see equation (2)) as a function of u {u = hu/kT). 
The physical conditions {log{T) = 6.3 and log{R) = —1.5) correspond to point 5 in Table 
1 and are close to the conditions at the base of the convection zone of the Sun. These 
monochromatic opacities have been obtained by interpolating OP mesh data (OPS data, 
Seaton private comunication) in Ne (electron density corresponding to logR — —1.5). Using 
these interpolated opacities to calculate kr, 7 and Qrad generate an error smaller than 2%. 
We estimated this error by comparing kr, 7 and Qrad obtained with the interpolated opacities 
to the values derived from the interpolation of 4 values of kr (or 7 or grad) calculated on the 
grid point. 

Prom the top and middle panels, one can see that new OP values are higher than OPl 
for both the individual and total monochromatic opacities. However increase in the specific 
elemental opacities is small for C, large for Fe and intermediate for the total opacities. This 
arises naturally from the changes in atomic physics .The inner-shell transitions included in 
the new OP data enhanced significantly the Fe monochromatic opacities, but lead to a small 
effect for C as expected (C is almost fully ionized). By definition, 7 is governed be the 
ratio of the individual to the total monochromatic opacities weighted by f,y. The weighting 
function f^, decreases rapidly at low and high frequency, damping all differences for these 
regions. 

At lower temperature, the importance of inner-shell decreases and for the physical con- 
ditions typical of envelopes the two sets of data are in good agreement (within 20%). The 
previous data were meant to be used for this purpose (Seaton et al. 1994). 

In the case of carbon the differences between the 2 ratios (NEW and OPl) will be 
reduced when we compare the 7 factors. This can be seen in the bottom left panel: k.^^ < 
and Kgf « Kgf"^ but < ^^^^^ 

1/ TOT TOT Kj^j Ktot 

For carbon, the increases in the total monochromatic opacities reduces significantly. As 
a consequence, the difference in grad cx /^i? x 7 will be smaller than in k.r. The changes in 7 
are partially compensated by those in kr 

In the case of iron, the fractional changes in the Fe monochromatic opacities are greater 
than those for Ktot, which makes 7 larger. This makes grad significantly larger compared to 
the OPl values. 

We calculated the acceleration for 15 elements present in the mixture for a range of 
{p — T) characteristic of the physical conditions found in different stellar models. The details 
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of theses models are given later in this section. For all the different structures, kr differs on 
average by less than 30% (see Table 3) with particular points differing by up to 38% (see 
Figure 1). For the acceleration, the results depend on the element and we can divide them 
in 2 groups. 

3.1.1. Light elements: C to Al 

The hghter metals follow the pattern seen in the carbon data. The differences in the 
acceleration are less than 40% with a rms smaller than 40% in a sense that the new data 
arc larger. These differences are dominated by the Rosseland mean and consequently the 
global effect is not. The monochromatic opacities for these elements are not significantly 
modified by the inclusion of the inner-shell transitions. The total monochromatic opacities 
are affected by the other abundant species for which the inner-shell transitions are important. 
As a consequence, the ratio of the two presents a difference that is balanced by the differences 
in the Rosseland mean opacity. The final results are closer than the Rosseland mean opacities 
themselves. 



3.1.2. Heavy Species: Si to Ni 

For the heavy elements, the accelerations increase by up to 80% with an rms between 
10% to 65%. These elements have 2 regimes. At lower temperatures, the increase in the 
monochromatic opacities is of the same order than in the total monochromatic opacities and 
the ratio of the two, which constitute the integrand in 7, stays relatively constant. Then the 
differences in the accelerations follow the trend of kr which does not differ a lot between 
the 2 datasets. At higher temperatures, where the contribution of the inner-shell transitions 
are very important, the ratio is not similar and 7 differs significantly {'^new > lopi)- The 
increase in the elemental monochromatic opacities is much more important than the rise in 
the total monochromatic opacities, kr is also larger, resulting a significant increase in the 
acceleration. 



3.2. OP and OPAL: Direct comparison for C. 

The monochromatic opacities for carbon obtained from the OPAL group allow us to 
directly gauge the effect of the difference in atomic data in calculating the radiative acceler- 
ations. 
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The OPAL data have been re-sampled and interpolated in order to match the frequency 
points used in the OP data. This does not introduce any significant error (less than 2% in 
all parameter). We estimated this error by recomputing the Rosseland mean opacity and 
comparing it to the value provided with the data by the OPAL group. In Table 4, the values 
in parenthesis are the recalculated values. We also reproduced Fig. 3 form Iglesias & Rogers 
(1995). We are confident that any error introduced by this procedure are negligible and 
allow a fair comparison. 

In Table 4 we compare the OP and OPAL Rosseland mean opacities values. The results 
for the Rosseland mean show a difference that docs not exceed 12%, in the sense that the 
new OP opacities are lower than OPAL (the differences between NEW and OPl range from 
6% to 30%). A detailed discusion on the comparison of Rosseland mean can be found in 
Badnell et al. (2004). In particular, the difference for the solar case is only 2.5% 

The results for 7 and grad are presented in Tables 5 to 8. The 7(C) factors differ by less 
than 12% when compared to OPAL and by less than 15% for OPl. The momentum transfert 
to the electrons (mte) does not change the results. The net results on the accelerations range 
from 5% to 12% when compared to OPAL without the momentum transfer to the electrons 
(mte) correction and from 2% to 17% when this effect is included. When compared to OPl 
the differences are less than 38% with mte corrections. 



3.3. OP vs OPAL: Other elements. 

The comparison with OPAL for other element can only be infer indirectly at the moment. 
The OPAL atomic data are not available. But from Seaton and Badnell (2004), their figures 
5,6,7,10 gives an insight of the expected differences. For these simple mixtures, one probes 
the difference of the element contribution in either dataset. As we have shown above, the 
difference on the carbon monochromatic opacities are small (Seaton and Badnell (2004) 
their figure 5). The same results is expected for O (Seaton and Badnell (2004) Fig. 6). The 
differences in S increase ( Seaton and Badnell (2004) Fig. 7) and become significant for Fe 
(Seaton and Badnell (2004) Fig. 10). All these results are for the condition of the 5 points 
from Table 1. An interesting point is that the new data produces larger monochromatic 
opacities for C, O and S compared to OPAL but smaller values for Fe. To extend the 
comparison, we have calculated the accelerations for different element for 4 stellar models 
and compare them to previous published works. 
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3.3.1. Models. 

For the other elements the four types of models used here are: 

1) The Sun at 4.57 Gyr; 

2) A model (R98) with Te// = lOOOOif and log{R) = -3 where R is defined (in previous 
OPAL works) as R — ^ with Tq — ^ where p is the mass density in g cm~^ and T the 
temperature in Kelvin. 

3) A M = 1.3M0 model at 70 Myr {T^ff = 6500X); 

4) A M = 1.5M0 model at 30 Myr (T^// = 7080/^); 

We picked regimes where the diffusion effects are known to be important. While levita- 
tion is a small effect in the Sun, we have precise data to compare theory and observations. 
The solar models therefore provide a useful point of comparison with other investigators. 
The other models are in regimes where the radiative levitation is most likely to play an 
important role; the different cases span a large domain of the (p, T) space. The second case 
is designed to mimic the physical conditions appropriate for hot horizontal branch stars or 
intermediate main-sequence stars of mass around 2.5 M©. The third and fourth cases are 
models of typical F stars where levitation is producing overabundances of Fe-peak elements 
for the slow-rotating Fm stars (T98a) . 

The results from T98a,T98b and R98 are taken from the articles directly, using a digital- 
ization of Figures 11 & 12 from T98a, Figure 1 from T98b and Figures 1 & 7 from R98 ( the 
original data for the structure, acceleration and monochromatic opacities were not available 
to us). 

The results for the Rosseland mean opacities are presented in Fig. 1. There are 4 panels 
representing the percentage difference in kr of the 4 stellar models studied. 

3.3.2. Results. 

The OPAL k,r have been calculated using the tables and the routines provided on 
the OPAL website . The results for the accelerations for the solar model at 4.57 Gyr 
are presented in Fig. 3. We have plotted the ratios using OP and OPAL data for 
several elements, where g is the gravity. There is a maximum difference of 10% in kr, but 



^ http : / / www-phys . Unl . gov/Research/ O PAL / 
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the difference in Qrad rises to 60% for Fe and up to 200% for C between the accelerations 
calculated with the new OP data and the data extracted from T98a, in the sense that the OP 
values arc smaller. For physical conditions appropriate for the base of the solar convection 
zone, the difference (defined as ^^^-^^^ ) is 20% for Fe and 150% for C. 

The accelerations for the second model with T^ff — 10000 K Log{R) — —3 are presented 
in Fig. 4. The Rosseland means agree within 15% and the radiative accelerations within 
50%. 

Figures 5 & 6 give the results for the next two models (M = 1.3 Mq and M — 1.5 Mq). 
The difference in kr (Fig. 2) reaches a maximum of 35% while the accelerations differ by 
less than 70%. 

For all cases, regardless of the type of star or the region of the (p, T) plane explored, 
the differences in the Rosseland mean opacities between OP and OPAL data are much 
smaller than the disagreements between the corresponding accelerations. As described in 
section 3.1.1, k,r cannot alone be held responsible for the large gap between g'^J^ and g'^^d^^- 
The combination of the discrepancies in kr and 7 defines the differences in the radiative 
acceleration. The accelerations in the different stellar model show large discrepancies for 
most elements. The largest ones are found in the solar model. 



3.3.3. General trends. 

In section 3.1 wc anticipated the differences in the acceleration for the elements heavier 
than carbon to be larger than for carbon itself. For the last three cases (T = lOOOOii', 
M = I.3M0 and M = I.5M0) the results seem to follow roughly the trend (sec Fig. 4,5 and 
6). However the differences are well within 50% with a rms of the difference smaller than 
30%. 

However for the solar case the differences in carbon accelerations are the largest. The 
difference for carbon range between few % to 250% through the solar structure while all of 
the other elements have differences within 60% (the regions near the center, small r/R^, are 
uncertain due to the difficulties to digitize this region). It is difficult to understand these 
results. Our direct comparison for carbon yielded a difference of 2% for case 5 where the 
the physical conditions are similar to those of the base of the convection zone of the solar 
model. The results from T98a give a difference of 150% at similar physical conditions. In 
other words, our accelerations using OPAL data are smaller than those presented by T98a 
for similar conditions. The points are not exactly at the same p — T, but it seems unlikely 
that this can explain such a large difference because the changes of T and p are very small. 
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We will address this issue in the next section. 



4. Discussion. 

4.1. Stellar Model differences. 

The first source of differences we can think of is the stellar code itself. We digitized 
figures from the published papers and attempted to construct equivalent physical cases for 
our models. Our models were derived from a different stellar evolution code than the one 
used by T98a, T98b and R98. In principle, differences in opacities and accelerations could 
therefore simply arise from differences in the structure of the models themselves (p, T, Xi). 
The total flux could also change the absolute scaling of the accelerations. The digitization 
process introduces errors. Even if the OPAL data were provided by the same authors (Carlos 
Iglesias), it could be possible that our OPAL data differ from theirs. Finally it is also 
possible that there are genuine differences in the calculated Qrad that arise from errors in 
either calculation. 

Modern solar models have very similar thermal structures for similar input physics, so 
we believe that structural differences are not an explanation for the differences in the solar 
case. When comparing our models with T98a (model H, Table 6), the base of the convection 
zone is at similar depth ( ^^(MoJ//, r98a) = O-'^l^G, = 0.7146), the temperature 

are in very good agreement {logiT^'^"') = 6.3343, log(T^'^'^) = 6.3360), the central density 
and temperature are within 1% and the composition is within few percent at the surface. 
We also checked that a 10% difference in the composition does not produce any significant 
difference in the acceleration. The constraints on solar models are the strongest compare to 
other stellar models. The two models use similar input physics (The Krishna-Swami T — r 
relation is used for the atmosphere, the energy generation routines are the same and both 
use the mixing length theory) . The equations of state differ but the effect is expected to be, 
and is, very small (for the center and the surface convection zone, p and T are within 1% in 
both models). Indeed for the interior solar conditions, the gas can be assimilated to a fully 
ionized ideal gas. It is very unlikely to have large difference in the structure. We calculated 
the partial derivatives of the acceleration at few points to gauge the uncertainties due to the 
structure. The partial derivatives are calculated as simple differences ^l^^g'^^x') ~ ^Aiog{x') 
where is T or p or r/i?*. Two of the variables are kept constant while calculating the 
partial derivative. 

At the base of the surface CZ of our solar model we have : ^^|^^^|p,r = ~5, ^^^^^^|r,r = 
0-5 and ^^f^|^UT = -1.22. 
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This implies that a deviation of 20% in the temperature or a factor of 7 in density or 
a difference of 112% in the radius or flux is needed in the model in order to invoke the 
structure as the source of differences in the radiative acceleration. As we mentioned above 
the differences arc less than 1% for the temperature, 2% for the central density and less than 
1% for the position of the base of the convection zone. 

For the second model the structure is just a vector defined identically in our model 
and the published one. It consists of a set of temperature points {log{T) =4.3 to 7.3) at a 
constant R {log{R) = —3, R = Tq ~ ■^). 

For the 2 other stars, the final model Tg// are identical to those from T98b for a similar 
age. So we expect the structures to be very similar; the effective temperature sets the 
convection zone depth, so the thermal structures should be similar at a fixed T^ff. 

Because the structure for the 2 last models are not as well constrained as the solar model, 
and because the differences in the acceleration are not too large it is a potential source of 
the differences between our calculation and the results from the literature. However, for the 
solar case, it is ruled out. 

The composition is also not the right candidate because they are equivalent to a solar 
composition in our model and in the 3 other models (2 last models are taken at the zero 
age main-sequence and by definition for case 2). For the solar case the composition does not 
change by more than 10% and we checked that this does not affect the acceleration. 

The errors due to the digitization are at the level of 3%, which is measurable but much 
less than the differences that we are seeing. 

The OPAL monochromatic opacities may be different from those used by the previous 
investigators. The data we received (Iglesias 2004, private communication) were re-created 
for 5 specific p — T conditions (Table 1). Despite the fact that we were able to reproduce 
Iglesias & Forest Figure 3a (Iglesias & Forest 1995) it is possible that the OPAL data used by 
Turcotte and others differ from ours. In the second case, where the structures are identical 
by definition, we found differences in the OPAL Rosseland mean opacities. The extracted 
value from the plot (R98 Figure 1) is hr = 3.45 at log{T) = 6.0 and log{R) = —3 while our 
value is = 3.72. 

Finally, we should not exclude the possibility that the two results are really different for 
the same physical conditions which could indicate an error in one or both calculations. 

We have been through the different possible sources of errors and we will not speculate 
further on these results given the number of inconsistencies. A direct comparison with the 
2 sets of data would be more reliable and we are expecting to be able to run it for other 
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elements than carbon and for a large portion of the p — T plane. We will be glad to provide 
our data to allow others to run parallel comparisons. 

4.2. Further discussion. 

The frequency resolution used to sample the monochromatic opacities, as already men- 
tioned by Seaton (1997,1999), is another important source of possible difference in the ac- 
celerations. The calculations of kr and 7 use a samphng technique (Seaton 1997; Leblanc, 
Michaud & Richer 2000) . All the cross sections are sampled at a specific frequency mesh and 
then integrated in order to calculate the Rosseland mean and 7. While a modest frequency 
resolution might be sufficient for the determination of /t^, higher resolution is required for 
the acceleration. Indeed, as discussed in Seaton (1997), 7 is much more sensitive to the 
frequency mesh than the Rosseland mean. The specific opacity will use a smaller number 
of line than kr. If the frequency mesh is bigger than some line width, the sampling might 
miss some lines. On average, it will have a small impact on the total opacity given the large 
number of lines already included. However it will be more crucial for the specific opacity 
which has less fines than the total opacity. 

While both OP and OPAL agreed on the convergence of kr within 2% (Rogers & Iglesias 
1992; Seaton et al. 1994) using 10^ points over the range of significant frequency, using an 
constant mesh (u), Seaton (1997) compared the effect of varying the frequency resolution 
from R = 10^ to R = 10^ to R = 10^ on 7. He showed that it could lead to a factor of few 
difference in 7, and varies depending on the element and the physical conditions. Seaton also 
addressed this issue in appendix B of his 1999 article. While the frequency resolution used 
in OPAL data (constant u mesh with 10^ points over the all frequency range) is adequate to 
determine with good accuracy the Rosseland mean opacity, it certainly will limit the accuracy 
of the acceleration when the composition changes during the evolution calculations. The new 
OP data use an equally spaced mesh (v) which samples preferentially the region where fi, is 
not too small (see Badnell et al. 2004) . This increase significantly the accuracy of 7 without 
requiring a prohibitive number of points. We have tested the results of the accelerations 
using 10^ points in v-mesh and compared with the acceleration using 10^ points in u-mesh. 
The results are in good agreement. We then expect more accurate accelerations derived from 
OP data than from OPAL. 
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5. Conclusion. 

The phenomenon of radiative levitation can induce significant changes in the surface 
abundances of stars, and the underlying theory that predicts levitation of some elements 
is based upon well-understood physics. However, quantitative theoretical predictions rely 
upon subtle details in the radiative opacities, and it has been difficult to estimate the error 
induced by uncertainties in the atomic physics and equation of state. The availability of 
two independent theoretical sources for opacity data (OP and OPAL) therefore provides an 
important test of systematic errors and their potential impact on this problem, as well as 
other aspects of the theory of stellar structure and evolution. 

We have calculated the radiative accelerations using the up-dated OP opacities which 
include all contributing inner-shell processes. The impact on grad is important and generates 
an increase by up to 80% compared to the data without the inner-shell. The light element 
are less affected than heavy species. The former changes by up to 40% while the later are 
affected to a 80% level. 

It is important to note that the difference in accelerations was significantly larger than 
the difference in the Rosseland mean opacities; we have traced through the interplay of 
factors that is responsible for the increased sensitivity of radiative acceleration computations 
to changes in the underlying atomic physics relative to the mean opacity. 

In the case of carbon, we were able to directly compare the OP and OPAL data for 
some cases. We found that the accelerations derived from the new OP data differ by less 
than 17 % compare to the accelerations derived from OPAL data. As expected, we found 
larger differences for heavier species in indirect comparisons with published literature values, 
ranging from 10% to 70% than the OPAL ones, depending upon the element and the physical 
conditions. 

However, these indirect comparisons also gave larger differences for C than the ones that 
we computed directly. In other words, our reconstruction of the C accelerations from the 
OPAL data differ from the published accelerations calculated using OPAL data in the sense 
that our OPAL accelerations are smaller. For the case that can be compared most directly 
(our solar case) , the differences were 2% while the differences with the published work reach 
150%. 

There are a variety of possible explanations. We had to digitize plots published for 
sample physical cases, and either digitization errors or differences in the published thermal 
structure relative to our models could contribute to the discrepancies. There could also be 
systematic differences in the accelerations for the same physical conditions that should be 
explored. In our view it would be most fruitful to simply perform a direct comparison rather 
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than to speculate. Until such a direct comparison can be made, it is reasonable to infer 
that different sources for opacities can produce differences at the factor of 2 to 5 level in 
the accelerations. We do note that the differences found for other elements were larger than 
the discrepancies between the two sets of carbon accelerations, suggesting that the relative 
trends (bigger differences for heavier species) are correct. 

The new OP data are in better agreement with the OPAL data than was the case for 
the older OP data (Badnell & Seaton 2003; Seaton & Badnell 2004; Badnell et al. 2004), 
and by extension we expect that the accelerations will also be closer as seen in the case 
of carbon. However, there arc real differences in the relative opacities of different species 
between the OPAL data and the new OP data, and there will be significant differences in 
the accelerations that we will discuss when we will have access to the OPAL data. 

There are two classes of direct astrophysical implications, one for the interpretation of 
stellar surface abundance anomalies and the other for tests of the opacities themselves. 

The differences in radiative accelerations will have important consequences on the micro- 
diffusion processes. The balance between thermal diffusion, gravitational settling and radia- 
tive Icvitation will be significantly modified. Lower values of Qrad could reduce the surface 
abundance anomalies, making it easier to reconcile the theoretical values with observations. 
It could also change the abundance profiles. Indeed the depth at which one species is sup- 
ported is directly linked to the balance between the 3 components of the diffusion processes. 
As a consequence, the Fe convection zone predicted in some F star models could lay deeper 
or even disappear. Levitation is only a perturbation in the solar case, and differences are 
unlikely to prodiicc significant modification in the calibrated solar-like models. However, as 
the effective temperature increases, the effects will increase too. Therefore, theoretical errors 
need to be accounted for when comparing observations and theory. 

More broadly, opacities play a crucial role in the theory of stellar structure and evo- 
lution. It has been difficult to establish the uncertainty in opacity calculations because of 
their complexity and the difficulty in obtaining direct measurements of opacity for plasmas 
under stellar interiors conditions of temperature and density. The sensitivity of radiative 
acceleration calculations to the monochromatic opacities, however, may provide a useful test 
of the opacities themselves in cases where the physical situation is relatively simple. By 
extension, such comparisons will be useful for establishing the quality of opacity data across 
the entire stellar evolution field. 

We would like to thank Carlos Iglesias for providing us with the carbon and total 
monochromatic opacity data from the OPAL project. We would also fike to thank Michael 
Seaton, Carlos Iglesias, and Anil Pradhan for useful discussions. FD thanks Georges Alecian 
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Part of this work was performed during stays by PD at the Observatoire de Paris, Meudon 
(Prance). PD is indebted to the INSU (CNRS, Prance) for support during those periods. 
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Fig. 1. — Percentage difference in k^i between OPAL and OP: upper left: Sun at 4.57Gyr, 
upper right: T^jf = lOOOOiT Log{R) = -3, lower left: M = l.SM© T = 6500K, lower right: 
M = 1.5M„ T = 7070ir 
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Fig. 2. — Monochromatic opacities at log(T)=6.3 and log(R)=-1.5: Left : for Carbon, Right 
: for Iron. Top : monochromatic opacities, Middle: total and Bottom: ^^i^. Dashed line 
-> OPl and solid line New OP 



- 23 - 



Solar model at 4.57 Gyrs 
— I — I 100 





0.5 
r/R. 



0.5 
r/R. 




- 'Fe' ' ' 


1 1 1 1 1 1 1 1 1 1 1 , 1 1 








< < < < 1 < < < < 1 < . < < 1 ' 



0.5 
r/R. 



Fig. 3. — Percentage difference in acceleration for tlie solar model. 
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Fig. 5. — Percentage difference in acceleration for M = I.SMq. 
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Fig. 6. — Percentage difference in acceleration for M = 1.5Mq. 
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Case 


Log{T) 


Log{R) 


Log{p) 


1 


6.0 


-1.5 


-1.5 


2 


6.0 


-2.0 


-2.0 


3 


6.0 


-2.5 


-2.5 


4 


6.0 


-3.0 


-3.0 


5 


6.3 


-1.5 


-0.6 



Table 1: (p-T) points at which individual monochromatic opacity and total opacities have 
been compared. For comparison, at the base of the solar convection zone in our model, 

Log{T) = 6.33 and Log{p) = -0.735. 



Element 


Z 


Number fraction 


H 


1 


9.071E-01 


He 


2 


9.135E-02 


C 


6 


3.770E-04 


N 


7 


9.913E-05 





8 


7.877E-04 


Ne 


10 


1.277E-04 


Na 


11 


2.271E-06 


Mg 


12 


4.039E-05 


Al 


13 


3.135E-06 


Si 


14 


3.769E-05 


S 


16 


1.722E-05 


Ar 


18 


3.518E-06 


Ca 


20 


2.434E-06 


Cr 


24 


5.047E-07 


Mn 


25 


2.608E-07 


Fe 


26 


3.359E-05 


Ni 


28 


1.985E-06 



Table 2: Compositon 
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Model 


(New - OPAL) 


(New - OPl) 


( OPAL - OPl) 


Sun 


1.3% 


29% 


29% 


10000 


6.2% 


6.5% 


9.5% 


M = I.SMq 


1.9% 


27% 


26.9% 


M = 1.5Mq 


4.8% 


17.4% 


17.6% 



Table 3: rms of the percentage difference for all points of each model, rms — y — ^ — 



log(T) 


log{p) 


^Present 


^OPAL 




(NEW-OPAL) 


(NEW-OPl) 


6.0 


-1.5 


58.75 


60.38 (60.37) 


44.63 


-2.7% 


24.1% 


6.0 


-2.0 


27.54 


28.83 (28.90) 


22.10 


-4.7% 


19.8% 


6.0 


-2.5 


10.29 


11.42 (11.52) 


8.93 


-10.9% 


13.2% 


6.0 


-3.0 


3.41 


3.65 ( 3.72) 


3.22 


-7.0% 


5.6% 


6.3 


-0.6 


34.89 


35.79 (35.76) 


24.61 


-2.6% 


29.5% 



Table 4: OPAL and OP kr and 7c. The percentage difference correspond to ^ opal — . 

The OPAL Rosseland mean has been recalculated using a sampling/interpolation 
procedure in order to have the same frequency point. This procedure has an error of less 
than 1% except for point 4 {log{T) = 6.0 log{p) = —3.0) for which the difference between 
the recalculated value (in parenthesis) and the value given by OPAL is 1.8%. 



-29- 



logT 


log{p) 


^og{Teff) 


r/R, 




logi-fOPAL-^ 


571 
7 


6.0 


-1.5 


3.762 


0.7146 


2.668 


2.664 


0.9 % 


6.0 


-2.0 


3.762 


0.7146 


2.583 


2.579 


0.9 % 


6.0 


-2.5 


3.762 


0.7146 


2.491 


2.476 


3.4 % 


6.0 


-3.0 


3.762 


0.7146 


2.406 


2.407 


-0.2 % 


6.3 


-0.6 


3.762 


0.7146 


2.239 


2.188 


11.1 % 



Table 5: Comparison of 7(C) obtained with OPAL and OP data (New). The effects of 
momentum transfer to the electron are not taken into account. 
Note: ^ = ^-^^^^^^^^^^ 



7 



logT 


log{p) 


log{Tef}) 


r/R^ 




log{^OPAL^ 




571 
7 


572 
7 


573 
7 


6.0 


-1.5 


3.762 


0.7146 


2.621 


2.625 


2.585 


-0.9% 


8.0% 


-8.8% 


6.0 


-2.0 


3.762 


0.7146 


2.526 


2.535 


2.496 


-2.1% 


6.7% 


-8.6% 


6.0 


-2.5 


3.762 


0.7146 


2.420 


2.424 


2.396 


-0.9% 


5.4% 


-6.2% 


6.0 


-3.0 


3.762 


0.7146 


2.322 


2.343 


2.299 


-5.0% 


5.2% 


-9.6% 


6.3 


-0.6 


3.762 


0.7146 


2.089 


2.060 


2.020 


6.5% 


14.7% 


-8.8% 



Table 6: Comparison of 7(C) obtained with OPAL and OP data (New and 0P5). The 
effect of momentum tansfcr to the electron is taken onto account. 

1\T + . '^71 — 7^eu)_^OPAL_ _ yNew_^OP5 _ g^^ _ ^OP_^OPAL 

iNOte. — jNew ) jNew ) jUPAL 
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logT 


log{p) 




r/R, 




iogig^f/"^) 


Sgi 

9 


6.0 


-1.5 


3.762 


0.7146 


4.077 


4.097 


-4.7% 


6.0 


-2.0 


3.762 


0.7146 


3.664 


3.692 


-6.7% 


6.0 


-2.5 


3.762 


0.7146 


3.145 


3.190 


-10.9% 


6.0 


-3.0 


3.762 


0.7146 


2.580 


2.629 


-11.9% 


6.3 


-0.6 


3.762 


0.7146 


3.423 


3.393 


6.7% 



Table 7: Comparison of the C acceleration obtained with OPAL and OP data (New). The 
accelerations do not take into account the effect of momentum tansfer to the electron. 
When the correction is applied, the the differences increase. 

'-9°' 
rjNew 



Xn^ n^ew ^OPAL 

Note: % = ^ .^1, 



9 9' 



logT 


log{p) 


logiTeff) 


r/R^ 






lo9{9?J^') 


9 


&92 
9 


9 


6.0 


-1.5 


3.762 


0.7146 


4.030 


4.058 


3.887 


-7% 


28% 


-32% 


6.0 


-2.0 


3.762 


0.7146 


3.607 


3.648 


3.494 


-10% 


23% 


-30% 


6.0 


-2.5 


3.762 


0.7146 


3.074 


3.138 


3.000 


-16% 


16% 


-27% 


6.0 


-3.0 


3.762 


0.7146 


2.495 


2.565 


2.460 


-17% 


8% 


-21% 


6.3 


-0.6 


3.762 


0.7146 


3.273 


3.266 


3.064 


2% 


38% 


-37% 



Table 8: Same as table 7 but the momentum transfer to the electron is taken into account 
as well as the scatering has been removed from the carbon mocrochromatic opacities 

, ,T Ar,. „New -OPAL r-„ -New _„OPl l„„ r,OPl „OPAL 

using the OP data. Note: ^ = ; ^ = ' ^^JL ] ^ = ' gJi^Ai 



